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"a^': 

^ ■ 1. Introduction 

^~~* . The present paper introduces a method of analysis of one dimensional sys- 

Zj , tems. Its application to study of Poincare recurrence time flow is considered. 

O ! The approach suggested consists of reducing the research of a given system's 

\^ ' orbits X = {xi)'^i to analysis of alternations of the monotone increase and 

decrease of higher order absolute finite differences, taken from X. Through 

^ . some special representation of finite orbits, we associate with X some nu- 

•^^ ! merical sequence, is called the conjugate orbit. We formulate two conditions, 

Q I when the study of the X is reduced to analysis of the special orbits' asymp- 

O ' totical intersections with some base set from numerical interval (0, 1) of zero 

Q. ■ Lebesgue measure. Hence, the approach can be characterized as an "asymp- 

Q\ . totical" analogy of Poincare's classical "section" method. 

^ ', The method distinguishes two cases: the continuous, when some numeri- 

^ p^' cal quantities p can take arbitrary values from interval (0, 1), and the discrete 

case, when all of them belong to a finite set. The continuous was applied in 

Ref. 1 to study of neural activity. In the present work, we are more interested 

rS^ . in the discrete case, which has an additional specific. Thus, as soon as we 

. . ! have stated the method's applicability, the mechanism generating time series 

.!^ I A, is replaced by the return map 7^, generating the conjugate series. While 

/\t ' the actual exact description of this mechanism may remain unknown (e.g., 

c3 ■ for earthquake time series or neuron spike trains), the function TZ does not 

depend on it and has a simple analytic shape. 

This approach is applicable to the study of various time series, arising in 
modern applied science (e.g., medicine, astronomy). The innovation consists 
of a new way to measure the fast oscillating time series - this measure is the 
Hausdorff dimension of some thin spaces A, to which the conjugate series are 
attracted. As an example of such an application, the computational analysis 
of the sequences of fractional parts is presented. We show that the results 
obtained testify on the possibility of Cantorian structure of our usual time 
flow. In this context, we recall some well known (qualitative) conclusions 
on irregular and intricate time flow, as indicated in Bergson's philosophical 
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treatments [0] and in some works on quantum physics ^, Q. 

2. Description of the Method 

2.1 Finite differences and conjugate orbits 

We deal with one dimensional, deterministic or stochastic systems, generating 
numerical sequences X = (a;,)^^ (0 < Xj < 1). We impose no restrictions on 
the system, and it may also possess different inner states changing with time. 
Below, we set two conditions on a given orbit X, under which it should be 
regarded as chaotic. We introduce the notion of conjugate (with X) orbit i', 
in terms of which the irregular nature of X can be established and studied. 
Let us first give some special representation of finite orbit X^ = (xj)^^^, 
which reflects its higher order differential structure. For the sequence X^ and 
number 1 < s < fc — 1 we let 



Ar=x„ A^> = \A^>-Ar~"\ {^ = 1,2,..., k-s). 



It is not difficult to obtain 






^t'^ = /^M-1 + E(-l)'" A(^^ - „,in (E(-l)^" AW) (1) 



0<i<k-s 
p=l ~ ~ p=l 



where 



. n A^''^ > A(*) ° 

^P^ = , aS' ; Ms) (it is supposed that ^ = 0) 



and 

/ifc,s = min{Aj^* : 1 < i < k — s} . 

From Eq. (1) it is easy to see, that if for a given m, < m < A; all of the 
finite binary sequences 

6[^6!{\...,6it (. = 0,l,...,m) 

as well as all of the numbers 

A^fc,!, /^fc,25 ■■■■,l^k,m ana iXi ,1^2 , . . . , ax^_^ 



are known, one can completely restore the initial finite orbit Xk = (xj)^^^. 
Hence, it follows from (1), we can consider that finite orbits with length k 
are given namely in the next special representation 



^k — [ri ,^2 , . . . ,r„ ,/^fc,i5 


/Ufc,2, • 


••,/ifc,m;Pl, P2, ■■.,Pk-m) (2) 


where 


■4'i 


(1 < s < m) (3) 


are the rational numbers, and 






l^k,i, l^k,2, ■ ■ ■,l^k,ni and pi. 


P2, •• 


. ,Pfc-m (where pi = aJ"^) 



are some numbers from interval [0, 1]. Here m = ruk and ruk tend to oo as 
/c ^ oo. It is easy to see, that after applying the recurrent procedure (1) the 
sequence X^ is completely recovered by (k- Now, we can introduce the basic 
tool of the method - the notion of conjugate orbit: we say, that numerical 
sequence u = {i's)'^i is the conjugate orbit, associated with the given orbit 
X = {xi)°Zi, if fo^ each s > 1 the terms z/^ are defined as follows 

u, = lim rf ) (= 0.6['^6i'^6i'^ . . .) 

fc— >oo 

were r^.'^^ are the r-coordinates from (2). 

Let us introduce a thin set B and its subsets Bk, are the base sets, men- 
tioned in Sec. 1. We consider the numbers < x < 1 represented in the form 
of binary expansion, 

oo 

X = 0.5,5253... (=^2-'=5fc, 4 = 0,1); (4) 

k=l 

further, we also operate with the corresponding binary sequences x, 

X = {5l,52,53,. . .,5n,. . .) . 

For a given natural k > 2 we define the subsets Bk of numerical interval 
(0, 1): Bk is the set all of those real numbers x G (0, 1) for each of which 

ni+i -ni<k (i = 1,2,3, ...) 

where rii denote all the consecuitive positions where the changes of binary 
symbol from (4) occur, 5^+^ = 1 — 5„.. Let also B be the union of all Bk 



{k > 2). All of the sets Bk, 2 < k < oo, have zero Lebesgue measure 
(see Ref. 1) (hereafter, in order to reduce some formulations, we use also the 
notation B^o (= B)). 

Let us now suppose, that the orbit X is such, that each of the sequences 

X-^ = i5['\5!^^...,5i:\...) (s = 0,l,...) (5) 

where s is fixed, has bounded lengths of series with the same binary symbol. 
In other words, we assume that for increasing sequence of indeces n^ , des- 
ignating all of those positions in natural series, where the changes of binary 
symbol occur, 6 \l,-^ = 1 — 5 %-^, we have 



n^^,-n^^<Ks<oo (^ = 1,2,...). (6) 

Hence, taking both numbers k and k — m large enough (for determinity it 
can be chosen m is equal to entire part of k/2, m = [k/2]), we can consider 
each of the sequences (5) as the binary expansion of some number Ug G Bk, 

or, in other words, each of the sequences {r^''^)'^-^ of the form (3) converges 
(as k -^ oo) to a number Ug from Bk- Here, 2 < K < oo, and according to 
(6) it should be taken K = sup {Ks '■ s > 1}. 

Continuous Case. Now, we asume also, that the quantities 

rrife k—mi^ 



converge to zero. 



hm /ifc = . (7) 

fc— >oo 



In such a way, the limitations we impose on the time series X are the follow- 
ing: 

(Ci): for every k>l the sequence X^ G B 

(C2): the quantities pi^k and jii^k are such, that: pk = o(l) as 

/c ^ 00. 



Then, one can see [|T|, these rwo restrictions imply that the transformed 
sequence (k from (2) tends to space B°° {= B x B x . . .), 

\\(k-B^\\ = o{l) (k^oc) (8) 

where ||.|| is the usual metric in Euclidean spaces i?^, ||x|| = {Y.i=ix'iY^'^- 
Particularly, (8) implies that for every s-th {1 < s < m) coordinate of Cfc we 
have 

lim Ct^ = lim rf ^ = u, e B 

k^oo fc— >oo 

- namely this relation was meant when we mentioned the "asymptotical" 
intersections (of transformed into a special form orbits) with a thin space. 

Discrete Case. In this case we deal with one dimensional time series with 
upper bounded terms and with restricted measurement accuracy, known in 
advance. If it is some number of the form 10"™, m > 1, then multiplying all 
of the terms of our series by this quantity, we will obtain some time series 
X, consisting of natural numbers bounded by a pregiven N, 

X = (ni,n2,n3,...,nfc,...), n^G [0,A^] {k = 1,2,3, ...) . (9) 

In other words, we deal with the random sequences ; particularly, such kind 
of sequences arise as a result of the reahzations of hazard games. However, 
we consider only a special class of sequences (9) - as above, we impose two 
limitations (below d > 1 is some natural number; we may assume d = 1) on 
the differential structure of the series X: 

(Di): for every A; > the sequence X'' G Bk 

{D2): for all large enough k > 1 and all i > 1 the quantities 
pf G{0,rf}. 

It should be noticed, that condition (Di) (and (Ci)) implies the dimension- 
ality limitations (see next section), while the (-D2), an analogy of (C2), is a 
requirement on certain regularity of the process, generating X. 

The restriction (-D2) means, that we consider such processes (9), which 
are reduced to binary processes X (with components Xi = 0,1). Let us 
call the binary sequence X, satisfying the conditions (Di) and {D2) as j3k- 
sequence and denote Bk the collection all of the numbers x G (0, 1) for 



which its binary expansion x is /3j^-sequence. From the definition, we have 
the next characteristic of these sets; below x„ = A^")(x) and A is the shift 
transformation (or the "tent map", see Ref. 6): 



Mx] 



2x < X < 1/2 

2a; - 1 1/2 < a; < 1 



Proposition 1 The number x G Boo iff for some numbers < Kp < cxo the 

relation 

|a;„-A(p)(x„)|>2-^^ (10) 

holds for every n > and p > 1. The number x G Bk, 2 < K < oo, iff' the 
inequality (10) holds with Kp = K and for every n > and p > 1. 

According to Poincare recurrence theorem (see, e.g., Ref. 6), the relation 

liminfp^ook- A(P)(a;)| =0 

holds for each < a; < 1 except some set of zero Lebesgue measure. In the 
case Kp = K the condition (10) means recurrence break (by analogy with 
"ergodicity break" from Ref. 7) for all points x„. Hence, all the sets Bk are 
found within exceptional set, mentioned in just quoted Poincare theorem. 
Note also certain resemblance of the relation (10) with the basic inequalities 
in Li - Yorke known theorems ^j. 

2.2 Conjugate attractor and return map 

Let us now denote A the closure of conjugate orbit {i's)%i , i-e. A C [0, 1] is 
the union of the set {i^s : s > 1} with the collection of all its cluster points. 
Then this orbit should be considered as a chaotic one, whenever A fills in 
either an interval or a Cantor set. We note, that for the most actual systems 
the constants from Eq. (6) for all s > 1 are upper bounded by the same 
number K >2 and, consequently, 

A C Bk and dim{A) < dimlBx) ■ (H) 

Hence, just mentioned criterion of chaosity simply means that ^ is a Cantor 
set. If A appears to be the same set for almost each orbit X, then we can 
conclude that A is the attractor of the system. The numerical values of 
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Hausdorff dimension of the sets Bk {2 < k < oo) can be obtained from the 
formula [Q] 

k 

z^^2'^im{Sk)^ = 1 where Zk{s) = J2 «"" (* > 0) 

n=l 

- hence, (11) imphes an upper estimate of Hausdorff dimension of attractor 
A. 

If we consider the numbers x from unit interval are given in the form (4), 
then it is not difficult to obtain the next statement (below, a(Bb = \a — b\ is 
the logical sum of the binary symbols a,b E {0, 1}): 

Proposition 2 If X is a binary sequence, then the return map 7?. : z/^ — *■ 
Uk+i of the conjugate orbit v has the form 

oo 

7^(a;) = 5]2-"(e„©e„+l) 

71=1 

where e„ are the coefficients of binary expansion of number x. 

One can find (in the operator form) the function TZ in Sinai's early construc- 
tions on weak isomorphism of dynamical systems [1^, |ll[]. For more details 
on this function see Ref. 9, where a fractal dynamical system is introduced. 

3. Sequences of Fractional Parts and Time Flow 

Some works in quantum mechanics 0] introduce the concept "Cantorian 
Space-Time" . The computational study of the fiow of our usual (but in some 
sense discretized) time, implemented in this section, confirms again the pos- 
sibility of such a concept. We are based on the method introduced - an 
approach completely different from those considered in other works on time 
fiow0,|T2|,|Il. 

Let us start from the statement and computational analysis of the follow- 
ing number-theoretical problem. Let us have two numbers < a, A < 1, a 
is irrational and let 

{«}, {2a}, {3a},..., {na},... (12) 

be the sequence of the fractional parts, generated by a (hereafter {x} denotes 
the fractional part of number x). We are interested in the behavior of the 
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time series 

X = (rifc)fcli where Uk = m^+i - m^ (^ > 1) ; 

here mk are all of the natural numbers, arranged in increasing order, rrik+i > 
rrik, such that 

{amfc}e(0,A). (13) 

According to H. Weyl theorem [^, |T^ the sequence (12) is uniformly dis- 
tributed in the interval (0,1), so that we have infinite number of m^, sat- 
isfying (13). If we take into attention the ergodicity of the transformation 

TT{x) = {x + a}, vr: (0,1)^(0,1) , (14) 

this statement can also be deduced (for a.e. irrational a) from Poincare's 
recurrence theorem. From numerous theoretical results on the sequences of 
fractional parts (see details in Ref. 14), we note the following Slater-Florec 
theorem: for given irrational a there exist two natural numbers p and q, 
such that for all fc > 1 we have Uk G {p,q,p + q}- We may also note 
a formula, obtained in Ref. 15 - it expresses the exact value of the total 
number of indeces ruk, satisfying (13) and not exceeding a given N, through 
the coefficients of some expansions of the numbers N and A. 

Our basic conclusion in the problem stated, deduced from computational 
experiments, can be formulated in the next form: for each irrational < a < 
1 there exist the numbers < A < 1 such that the conjugate to X orbit u 
is asymptotically close to a Cantor set; and conversely, for each < A < 1 
there exist irrational numbers < a < 1, such that the conjugate to X orbit 
z/ is asymptotically close to a Cantor set. 

Some of these computational results, on the Fig. 2 are presented. Here, 
the numerical value A = 0.13 and several values of parameter a are consid- 
ered. In dependence of control parameter a the conjugate orbits p demon- 
strate all possible types of motion: the periodic, when they are attracted to 
some finite set, Cantorian, when the attractor A has Cantorian self-similar 
structure, and the completely chaotic motion, when the orbits fill in an in- 
terval from (0,1). For instance, for a = 0.1250002 the Cantorian attractor 
in Fig. 2.b is shown. The graph Fig. 2. a has been obtained for the value 
a = 0.124999; in this case the orbit u consists of 13 straight lines (some of 
them are not shown). In some respect, the behavior (in dependence of con- 
trol parameter a) of conjugate orbits in this problem reminds Feigenbaum's 



transitions (see e.g., Ref. 6), and this needs further detailed computational 
study. 

To end, let us consider the following interpretation of the above intro- 
duced number-theoretical problem; in particular, this will explain the title 
of present section. Let us have in a two dimensional plain a circle E and a 
point S, situated outside of disc, bounded by E. Let u be the angle, under 
which the circle E from the point 5* is seen, and let also L be an arc of E 
with the length \L\ = X < u with the center at the point of intersection of 
circle E with the rectilinear segment, connecting its center with S. Let us 
assume that the point M G -E is rotated in discrete time T = {1, 2, 3, . . . .} 
with the frequency a. Then, the consecuitive moments of falling the point M 
in the arc L constistute a sequence of natural numbers nik G T, coinciding 
with the sequence X from Eq. (13). Also, it is not difficult to see, that the 
sequence of fractional parts (12) coincides with the sequence of iterates of 
the function 7r(x) from Eq. (14), i.e. 

{an} = 7r(")(a) {n > 1) 

and, consequently, the numbers mj^ represent the Poincare recurrence (of 
point M into the arc L) time. Now, let us interpret E as "Earth" and S as 
" Sun" . Then, replacing the continuous rotation of Earth by discrete motion, 
we obtain that rrik appears to be an approximate of the flow of our usual 
time. In such a way, the above described computations demonstrate, that 
in principle, it is possible that our time possesses the Cantorian structure. 
The numerical values of the quantities a, A, and u can be specified for this 
"astronomical" case, but here we do not dwell on this point. 

4. Some Remarks and Discussions 

Let us discuss two problems, relating to the brain activity. In the paper 
Ref. 1, through the method above, applied to experimental data of actual 
neurons, the Cantorian strucuture of brain activity has been deduced. The 
method permits further development |^ by means of introducing some dy- 
namical system JF = [A, T, ^) on the Cantor spaces A. The evolution opera- 
tor T is easly expressed through the return map TZ, while the attractor A and 
invariant (singular) measure /i remain the main quantities, distinguishing dif- 
ferent processes. In practical applications, the measure /i can be calculated 
through frequency analysis, accepted in mathematical-linguistic works (see, 

9 



e.g. Ref. 5). Hence, we can operate with the theoretical-information charac- 
teristics of system T . Whether there exists an analogy of known Billingsley- 
Eggleston formula |TD 



Dimension = Entropy, 

relating to Eggleston thin spaces, for the fractal dynamical systems JF? Note 
also Mandelbrot's work |jT6|], who considered the applications of this formula 



to problems of turbulence. 

We know [|I^, that there exists some mechanism of time perception in 
the brain. Its activity, it follows from our previous results |Q , must be found 
in some thin (multidimensional) Cantor space. How to relate the measur- 
able spaces J-'t and J-'p, which correspond to time flow and time perception, 
with each other? Particularly, can we compare the values of dini{At) and 
dim{Ap)7 

It seems, these problems can be resolved after detailed computational 
frequency analysis, in order to compute the invariant measures /i (and hence 
another characteristics, say, entropy) of both time flow and brain activity. 
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Figure Captions 

Figure 1. The graph of inverse function Ai{x) = TZ^^{x). The graph is 
constructed through 15,000 points, computed by means of random numbers 
generator. From definition of the function 7^ it is easy to see that A1 is a 
multivalent function, defined on the interval (0, 1). The graph of the return 
map TZ can be obtained as a result of symmetric reflection of this "fractal 
letter M" in respect of diagonal (0,0), (1, 1). 

Figure 2. Two different types of conjugate orbit's behavior in dependence of 
control parameter a. The graph a) demonstrates the case when the attractor 
A consists of the finite set, while in the graph b) these orbits are attracted to 
some Cantor set. The graph c) shows the self-similar structure of the Cantor 
set A from the previous graph. 
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